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INTRODUCTION, 


Transport  and  diffusion  of  dust  and  nonbuoyant  smokes  are  simulated  in  a 
realistic  manner  by  the  particle  model  originally  developed  by  Ohmstede  and 
Stenmark1  2  at  the  US  Army  Atmospheric  Sciences  Laboratory.  However,  some 
smoke  rounds  release  their  contents  explosively,  and  others  involve  a  burning 
substance,  such  as  white  phosphorus.  Chemical  agents  also  may  be  released  in 
such  a  manner  that  they  have  an  initial  buoyancy.  The  consequence  of  these 
processes  is  a  buoyant  cloud  that  is  taken  Into  account  by  a  new,  simple 
algorithm  included  in  the  particle  model.  For  the  models  discussed  in  this 
paper,  a  particle  represents  a  small  parcel  of  air  plus  obscurant  or  agent. 
Each  parcel  is  so  small  that  expansion  with  time  of  individual  parcels  is 
insignificant  relative  to  the  scale  of  the  turbulence. 

This  method  uses  the  difference  in  temperature  between  an  individual  particle 
and  the  ambient  atmosphere  to  generate  a  vertical  acceleration.  This  method 
takes  into  consideration  entrainment  at  the  top  and  sides  of  the  cloud,  along 
with  ingestion  of  unheated  material  at,  or  loss  of  heat  by  contact  with,  the 
ground  during  initial  release  and  during  subsequent  ground  contact.  A  series 
of  sample  computer  runs  suggest  that  the  particle  model  with  the  new  algorithm 
can  model  realistically  transport  and  diffusion  of  buoyant  material. 


BUOYANCY 

Buoyant  Velocity 

The  particle  model  for  passive  materials  uses  a  method  based  on  the  Langevin 
equation  to  compute  particle  displacements.  In  very  simplified  terms,  this 
technique  assumes  that  at  a  given  time  step  a  particle  moves  as  a  result  of  a 
"remembered"  portion  of  its  velocity  at  the  previous  time  step  plus  a  velocity 
arising  from  an  acceleration  due  to  a  random  force.  Gifford3  and  Ohmstede  and 
Stenmark2  provide  a  much  more  complete  description  of  the  passive  model. 
Briefly,  in  the  form  of  a  simple  equation. 


w'(t)  *  w'(t  -  1)  *  R  +  w~(t) 


(1) 


1U.  0.  Ohmstede  and  E.  B.  Stenmark,  1980,  ”A  Model  for  Characterizing 

Transport  and  Diffusion  of  Air  Pollution  in  the  Battlefield  Environment/ 
Proceedings  of  Second  Joint  Conference  on  Applications  of  Air  Pollution 
Meteorology  and  Second  Conference  nn  Industrial  Meteorology,  Amer  Met  Soc, 
Bost&Ti,  MA - - 

2W.  D.  Ohmstede  and  E.  B.  S*.  mark,  1981,  "Parameterization  of  the  Dispersion 
of  Battlefield  Obscurants,  Proceedings  of  the  Smoke/Obscurants  Symposium  V, 
Adelphl,  MD 

3 F.  A.  Gifford,  1981,  Horizontal  Diffusion  in  the  Atmosphere:  a  Lagranglan- 
Dynamical  Theory,  Los  Alamos  National  Laocratory,  Keport  LA-»bb/-MS,  UC-J4b, 
Los  Alamos,  M 
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where  w'(t)  s  velocity  at  time  t,  w'tt  -  1)  =  velocity  at  the  previous  time 
t  -  1,  w"  =  velocity  arising  from  a  random  acceleration  (w"  =  n  *  (1  - 
R2)1'2  *  for  a  unit  mass,  where  n  =  a  random  acceleration),  and  R  is  a 
Lagrangian  correlation  function  with  a  Lagrangian  time  scale  relating  the 
velocity  at  time  t  with  that  at  t  -  l.2  3 

Ue  account  for  buoyancy  by  adding  a  third  velocity  term  (w'~)  on  the  right  so 
that 


w'(t)  =  w'* { t  -  1)  *  R  +  w"(t)  +  w'"(t) 


(2) 


The  equation  for  w"'(t)  is  based  on  the  simple  expressions  found  in  many 
meteorology  texts  and  references,  such  as  Haltner  and  Martin,14  Priestly,5  and 
more  recently  Ludlam.6  The  expression  used  here  may  be  considered  as  a  simple 
parcel  equation  in  which  an  acceleration  is  caused  by  the  difference  in  den¬ 
sity  between  the  parcel  {or  particle)  and  the  ambient  atmosphere.  This  dif¬ 
ference  may  be  related  to  the  corresponding  temperature  difference,  assuming 
the  pressure  is  the  same  Inside  and  outside  the  parcel.  With  T  =  ambient 
temperature,  TD  =  parcel  temperature,  and  g  =  acceleration  of  gravity  we  have, 
letting  w'"'(t)  be  replaced  by  , 


dw'~  =  g  *  C(Tp  -  Ta)/T#]  *  dt 


(3) 


2W.  D.  Ohmstede  and  E.  B.  Stenmark,  1981,  "Parameterization  of  the  Dispersion 
of  Battlefield  Obscurants,"  Proceedings  of  the  Smoke/Obscurants  Symposium  V, 
Adel  phi,  M) 

3F.  A.  Gifford,  1981,  Horizontal  Diffusion  in  the  Atmosphere:  a  Lagrangian- 
Dynamical  Theory,  Los  Alamos  National  Laboratory,  Report  LA-«bb/-MS,  UC-34b, 
Los  Alamos,  NM 

•♦G.  J.  Haltner  ano  F.  L.  Martin,  1957,  Dynamical  and  Physical  Meteorology, 
McGraw-Hill  Book  Co.  Inc.,  New  York,  NY 

5C.  H.  B.  Priestly,  1959,  Turbulent  Transfer  in  the  Lower  Atmosphere, 
University  of  Chicago  Press,  Chicago,  IL 


6F.  H.  Ludlam,  1980,  Clouds  and  Storms,  the  Behavior  and  Effect  of  Water  in 
the  Atmosphere,  Penn  State  University  Press,  University  Park,  PA 


Integrating  from  tp  to  t. 


w'~ 


9  *  C(Tp 


-  V'V  *  11  -  V 


(4) 


For  this  very  preliminary  investigation  we  made  the  further  simplification  of 
setting  t  -  t0  equal  to  1  s,  giving 


w 


>  +  A 


9  *  <Tp  -  Ta)/Ta 


(5) 


To  account  for  the  "aerodynamic"  drag  of  the  particle  against  the  surrounding 
atmosphere,  we  have  introduced  a  drag  coefficient  (k),  which  we  assume  only 
acts  on  the  additional  velocity  arising  from  buoyancy,  giving 


=  g  *  C(Tn  -  T  1/TJ  *  (1  -  k)  . 

P  9  3  16) 


Since  a  rising  parcel  tends  to  conserve  its  potential  tenperature  (©),  to  be 
strictly  correct  we  should  replace  T_  with  0_.  However,  the  use  of  0  would 
require  the  computation  of  additional  intermediate  variables  (for  example, 
atmospheric  pressure  at  each  particle  height)  thereby  increasing  computation 
time.  An  essentially  equivalent  method  is  to  use  Tp  along  with  the  dry  adia¬ 
batic  lapse  rate  to  account  for  adiabatic  cooling.  Therefore,  in  this  pre¬ 
liminary  study  we  have  used  ordinary  tenperature,  decreasing  it  according  to 
the  dry  adiabatic  lapse  rate  (0.0098  K  m-1 )  for  increasing  altitude.  For  a 
giverv  time  step,  when  the  buoyancy  routine  is  entered,  the  initial  particle 
temperature  T_(t  -  1)  Is  cooled  as  a  result  of  the  rise  of  the  particle  during 
the  previous  time  Interval  [(t  -  1)  -  (t  -  2)].  This  gives 


T  (t)  =  T  (t  -  1)  -  0.0098  *  [Z(t  -  1)  -  Z(t  -  2)]  , 

P  P  (7) 


where  Tp(t)  Is  the  value  after  adiabatic  cooling  but  before  consideration  of 
other  processes  that  modify  particle  temperature,  Z  =  height,  and  t  -  1,  t  -  2 
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refer  to  the  previous  time  steps.  Note  that  this  Tp(t)  is  computed  at  the 
beginning  of  the  buoyancy  routine.  Thus,  Z(t  -  1)  represents  the  height  at 
this  point  in  the  routine.  We  want  to  calculate  the  cooling  that  occurred  as 
a  result  of  the  rise  of  the  particle  during  the  previous  time  step,  that  is, 
the  rise  to  Z ( t  -  1)  from  Z{t  -  2). 


Entrainment 

Entrainment  is  handled  in  a  conceptually  simple  way,  which  assumes  that  buoy¬ 
ancy  of  particles  in  a  central  “column"  decreases  proportionally  to  their 
distance  from  the  centerline.  Outside  this  central  column,  the  particles  are 
assumed  to  lose  most  of  their  buoyancy;  that  is,  AT  =  Tp  -  Ta  becomes  small. 
The  same  holds  for  the  top  few  particles  that  represent  the  thin  uppermost 
layer  of  the  cloud. 

At  a  given  time  step  the  cloud  of  particles  is  divided  into  layers  having  a 
thickness  (aZ)  of  2  m  from  the  surface  (Z  =  0)  to  10  m,  and  a  aZ  of  5  m  above 
10  m.  The  mean  values  of  X  and  Y  (Y,  Y)  and  the  standard  deviations  (ox, 

Oy)*  are  computed  for  each  layer,  except  that  ox  and  are  set  to  zero  if  <  3 

particles  are  found  in  a  given  layer.  X  and  T  are  set  to  zero  if  a  layer 

does  not  contain  any  particles.  Figure  1*  shows  a  simplified  sketch  of  such  a 
cloud  in  two  dimensions  (X,  Z).  For  clarity,  only  a  few  layers  are  shown,  all 
of  which  have  the  same  thickness.  In  the  figure  the  central  region  is  out¬ 
lined  by  a  thicker  line.  This  "buoyant"  region  has  a  boundary  defined  by  the 
surface,  +ox  for  each  layer,  and  the  bottom  cf  the  top  layer  that  contains  s  3 
particles.  In  the  computer  algorithm  the  horizontal  cross  section  of  the 
region  has  the  form  of  an  ellipse  defined  along  the  X  and  Y  axis  by  +ox  and 
icy 

Outside  this  inner  region,  aT  is  modified  by  a  preset  proportion  depending  on 
the  value  of  an  entrainment  coefficient,  kg  (>  0.8).  Using  ke,  the  temper¬ 
ature  of  a  particle  is  calculated  by  replacing  Tp  with  Tp  -  at  *  kg. 

Within  this  "buoyant"  region,  AT  is  retained  at  each  time  step,  depending  on 

the  location  of  each  particle  relative  to  Y,  Y  of  each  layer.  A  particle  at 
X,  Y  suffers  no  loss  of  temperature.  For  other  values  of  X,  Y  the  retention 
of  temperature  depends  on  the  distance  from  the  origin  (Y,  Y)  of  the  ellipse 
for  each  layer.  The  percentage  of  AT  kept  falls  off  approximately  as  a 
parabola,  with  a  minimum  value  of  AT  *  (1  -  ke)  at  +ox,  +<y  Figure  2  shows  a 


*ox  and  o  are  the  respective  sanple  standard  deviations  of  the  x  and  y 
components a of  the  particle  locations  for  each  layer  of  arbitrary  thickness 
aZ.  These  quantities  are  not  the  standard  deviations  found  in  many  Gaussian 
plume  models. 

figures  are  presented  at  the  end  of  text. 
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schematic  of  the  scheme  for  a  given  layer,  in  only  two  dimensions  for  clarity. 
The  equation  defining  the  buoyant  region  has  variables  X,  Y,  X,  Y,  ox,  oy, 
and  we  let  f  equal  a  confuted  proportion  factor  with  values  from  0  to  1.  This 
gives 


f  =  1  -  (X  -  X)2  /  o2  -  (7  -  Y)2  /  J 

*  y 


(8) 


where  the  difference  between  the  mean  and  the  variable  is  not  allowed  to 
exceed  the  relevant  standard  deviation  (therefore  0  <  f  <  1).  Note  that  f  is 
not  computed  if  a  =  o  -  0;  that  Is,  If  the  number  of  particles  in  a  layer  is 
<  3.  Using  f  ana  ke/we  calculate  the  temperature  of  a  particle  in  the  cen¬ 
tral  column. 


Tp  =  T  -  (1  -  p)  *  AT 


(9) 


where  T  is  the  previously  calculated  particle  temperature  for  the  current  time 
step,  and  p[=l-ke*(l-f)]1s  the  portion  of  the  temperature  difference 
(at)  retained  (see  figure  2). 

Ingestion  of  Nonbuoyant  Material  and  Contact  With  Ground 

An  explosive  release  or  a  strongly  buoyant  quasi -continuous  release  will  cause 
the  ingestion  of  some  surface  and,  for  an  explosive  release,  subsurface  mate¬ 
rial.  This  additional  dust  and  other  material  will  add  a  nonbuoyant  component 
to  the  cloud.  In  the  real  world  some  of  this  material  may  be  heated  by  con¬ 
tact  with  the  air  and  smoke  (or  agent)  mixture,  which  in  turn  would  be  cooled 
by  the  same  contact.  Also,  some  of  the  cloud  material  may  cool  as  a  result  of 
contact  with  the  surface  during  Initial  release.  The  model  simulates  both  of 
these  processes  by  removing  a  part  of,  or  all  of,  the  buoyant  energy  (that  is 
AT)  from  some  of  the  released  particles,  according  to  a  unitary  random  number 
generator. 

The  original  model  assumes  that  particles  that  hit  the  ground  are  reflected 
perfectly.  Stenmark  (private  communication)  had  suggested  that  for  a  buoyant 
particle  (AT  >  0),  contact  with  the  ground  after  release  could  be  modeled  by 
having  It  lose  all  of  Its  extra  heat  (AT  =  0).  The  new  algorithm  modifies 
this  idea  slightly,  so  that  less  than  all  of  its  extra  temperature  may  be 
lost: 


AT  *  k.  *  (T  -  T  ) 


where  kg  is  a  coefficient  of  "heat  loss,"  Tp  =  temperature  of  the  particle 
before  contact  with  the  surface,  and  T$  =  surface  temperature.  Note  that  if 
Tp  <  T$,  the  algorithm  would  add  heat  to  the  particle,  since  T  is  subtracted 
from  Tp  in  the  calculation  of  the  "final"  Tp  for  the  given  time  step.  This 
computation  of  <‘.T  from  ground  contact  occurs  after  the  calculation  of  w'"', 
thus  changing  Tp  for  the  next  time  step. 


Other  Cooling 

Aside  from  entrainment  and  interaction  with  the  ground,  it  is  assumed  that  the 
difference  in  temperature,  \T  =  Tp  -  Ta,  decreases  exponentially  with  time. 
If  we  let  c  be  a  cooling  coefficient,  we  replace  aT  by  aT  *  (1  -  c),  where  c 
is  an  adjustable  parameter  ranging  from  0  to  1.  In  the  buoyancy  routine  this 
exponential  decrease  occurs  after  the  computation  of  w"',  thus  modifying  the 
particle  temperature  for  the  next  time  step. 


Statistics 

A  separate  subroutine  calculates  the  statistics  'Y,  Y,  o  ,  and  ->  used  to 

x  y 

compute  the  effect  of  entrainment  on  the  particle  temperature.  These  statis¬ 
tics  are  found  after  values  of  X,  Y  are  computed  for  all  particles  released. 
Therefore,  the  statistics  for  a  given  time  step  are  those  computed  for  the 
previous  step.  For  example,  the  proportion  factor  f  for  time  t  is  calculated 
using  X,  Y,  etc.,  for  time  t  -  1. 


Ambient  Temperature  Profile 

The  temperature  of  the  ambient  atmosphere  Is  assumed  to  vary  with  height 
according  to  a  simple  two-part  lapse  rate.  For  heights  up  through  a  prese¬ 
lected  altitude,  the  lapse  is  dry  adiabatic;  for  greater  heights  the  routine 
uses  an  arbitrary  lapse.  Since  the  particles  generally  do  not  rise  above 
several  hundred  meters,  a  more  complex  temperature  profile  would  not  produce 
any  significant  change  in  ambient  temperature  for  a  given  particle  height. 


COMPUTER  ALGORITHMS  AND  ROUTINES 

All  Input  variables  and  parameters  are  entered  by  the  name  list  $PDATA,  and 
all  have  default  values.  Table  1*  lists  the  input  variables  and  parameters 
for  buoyancy,  along  with  brief  explanations  and  the  respective  default 
values.  It  is  anticipated  that  in  the  future  one  or  more  of  these  quantities 
may  be  computed  instead  of  being  entered  as  input  data. 


♦Tables  are  presented  at  the  end  of  the  text. 


Ohrastede  and  Stenmark1  2  describe  the  computational  methods  for  the  model  for 
passive  materials.  Here  we  only  discuss  the  additional  routines  that  have 
been  added  to  account  for  buoyant  materials. 

The  routine  for  computing  statistics  has  the  form  of  a  subroutine  (STAT) 
called  from  the  main  program  (DRIVER)  at  the  end  of  the  computations  for  each 
time  step.  The  computed  statistics,  used  in  the  following  time  step,  include 
mean  values  and  standard  deviations  for  the  two  horizontal  coordinates. 
Standard  formulas  and  procedures  are  used  and  will  not  be  discussed  here. 

The  temperature  of  the  ambient  atmosphere  (TAIR)  at  a  given  particle  height 
(Z)  is  calculated  by  a  simple  function  (AIRFCN)  that  uses  Z  and  surface  tem¬ 
perature  (TSURF)  as  input.  If  Z  is  less  than  or  equal  to  a  preset  height 
(HITE),  the  temperature  equals  TSURF  -  ADLAPS  *  Z,  where  ADLAPS  is  the  lapse 
rate  (normally  dry  adiabatic)  for  the  bottom  layer  of  atmosphere.  For  Z 
greater  than  HITE,  the  temperature  is  given  by 


TAIR  =  TSURF  -  ADLAPS  *  HITE  -  LAPSE  *  (Z  -  HITE)  , 

(ID 


where  LAPSE  is  the  lapse  rate  for  the  top  layer. 

The  main  buoyancy  routine  is  part  of  the  subroutine  (STEPIT)  for  computing 
individual  particle  displacements  for  each  time  step.  This  algorithm  added  to 
STEPIT  is  shown  in  figure  3  in  the  form  of  a  flow  chart.  A  written  explana¬ 
tion  would  be  unnecessarily  confusing  and  tedious. 


SAI«»LE  OUTPUT 

Output  from  the  model  in  tabular  and  graphical  form  is  presented  for  both 
instantaneous  and  quasi -continuous  releases,  with  and  without  buoyancy.  The 
burn  time  for  the  quasi -continuous  cases  Is  30  s,  and  in  all  cases  the  data 
are  shown  for  a  model  time  of  150  s.  The  tabular  data  in  tables  2  through  5 
provide  (1)  the  number  of  particles  in  each  of  30  layers  up  to  a  height  of 


1 W .  0.  Ohmstede  and  E.  B.  Stenmark,  1980,  "A  Model  for  Characterizing 
Transport  and  Diffusion  of  Air  Pollution  In  the  Battlefield  Environment," 
Proceedings  of  Second  Joint  Conference  on  Applications  of  Air  Pollution 
Meteorology  and  Second  conference  on  industrial  Meteorology,  Amer  Met  Soc, 
Boston,  MA 

2W.  D.  Ohmstede  and  E.  B.  Stenmark,  1981,  "Parameterization  of  the  Dispersion 
of  Battlefield  Obscurants,"  Proceedings  of  the  Smoke/Obscurants  Symposium  V, 
Adel  phi,  MD 
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135  m,  (2)  the  means  and  standard  deviations  for  layers  with  four  or  more 
particles,  and  (3)  the  height  of  the  highest  particle.  Figures  4  through  7, 
presenting  the  same  cases,  present  particle  plots  for  the  X,  Z  plane.  Figures 
3  and  9  show  Y,  Z  plots  for  the  Instantaneous  releases. 

Throughout  the  remainder  of  this  discussion,  reference  should  be  made  to  the 
appropriate  tables  or  figures.  The  maximum  spread  of  the  cloud  occurred  with 
an  instantaneous  release  with  buoyancy  (figure  4),  followed  by  the  buoyant 
quasi -continuous  case  (figure  6).  Without  buoyancy  the  clouds  for  the  instan¬ 
taneous  (figure  5)  and  quasi -continuous  (figure  7)  releases  are  fairly  simi¬ 
lar,  the  former  showing  a  little  more  dispersion.  The  respective  tables 
readily  illustrate  the  greater  vertical  extent  of  the  buoyant  clouds.  Instan¬ 
taneous  releases  (tables  2  and  3)  yield  more  extensive  clouds  than  their 
quasi-continuous  counterparts  (tables  4  and  5).  In  the  Instantaneous  cases, 
all  particles  have  a  history  of  150  s,  while  the  model  time  for  the  quasi- 
continuous  particles  ranges  from  120  to  150  s.  Therefore,  one  could  expect 
the  Instantaneous  cloud  to  grow  somewhat  larger  than  the  quasi-continuous  one, 
other  conditions  being  the  same. 

The  turbulence  parameters  on  release  are  assumed  to  be  those  of  the  ambient 
atmosphere,  as  for  a  nonbuoyant  release.  Any  "extra"  turbulence  arising  from 
a  heated  (buoyant)  release  is  assumed  to  be  accounted  for  by  the  buoyancy 
routines.  For  example,  the  net  radiation  Index  Is,  as  before,  set  at  0,  the 
mean  windspeed  at  5  m  is  2  ms-1,  the  power-law  exponent  for  windspeed  is  0.2, 
and  the  10-min  standard  deviations  of  the  wind  azimuth  and  elevation  angles 
are  10°  and  3.5°,  respectively. 

Figures  8  and  9  show  the  greater  crosswind  spread  of  the  buoyant  cloud  in 
spite  of  the  cloud  being  oriented  in  the  along-wfnd  direction  (mean  crosswind 
component  of  wind  =0).  In  the  real  world  this  increase  in  lateral  dispersion 
may  arise  from  Increased  turbulence  at  greater  heights  (mixing  lengths  are 
longer).  Buoyant  particles  rise  to  a  given  height  more  quickly,  thereby 
giving  the  larger  scale  turbulence  more  time  to  disperse  them.  In  the  model, 
dispersion  depends  on  the  product  of  the  variance  of  vertical  velocity  and  the 

Lagrangian  time  constant  (o^  TL)  .  With  the  given  input  for  these  sample 

runs,  changes  very  slowly  with  height,  but  T^  increases  relatively  rapidly, 

allowing  more  dispersion  of  particles  at  greater  heights.  Thus,  buoyant 
particles  disperse  at  a  faster  rate  and  have  more  time  at  higher  levels. 

Particles  that  contact  the  ground  lose  their  buoyancy  (AT  set  to  0),  and  a 
preset  percentage  of  particles  on  release  are  considered  to  be  ingested  non¬ 
buoyant  material,  or  to  have  lost  their  extra  energy  by  contact  with  the 
ground.  Therefore,  the  downwind  travel  distance  of  particles  very  near  the 
ground  should  be  similar  for  the  buoyant  and  nonbuoyant  run  for  each  type  of 
release.  Comparison  of  figure  4  with  5,  and  6  with  7,  shows  that  very  close 
to  the  ground  the  most  traveled  particles  (to  the  right)  occupy  similar 
positions  In  the  X  direction,  as  do  the  least  traveled  particles  (to  the 
left).  For  both  Instantaneous  and  quasi-continuous  releases,  the  fastest 
rising  particles  for  the  buoyant  runs  moved  further  downwind.  This  greater 
movement  was  expected,  since  windspeed  Increases  with  height  in  the  model. 
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PRELIMINARY  SENSITIVITY  ANALYSIS 

A  preliminary  analysis  was  undertaken  to  obtain  a  rough  estimate  of  the  sensi¬ 
tivity  of  the  particle  model  with  buoyancy  to  several  of  the  factors  listed  in 
table  1.  An  instantaneous  release  was  chosen,  since  it  produced  the  greatest 
reaction  to  the  addition  of  buoyancy.  Output  for  a  model  time  of  15C  s  was 
compared  with  the  results  shown  in  table  2.  Tables  6  through  14  present  the 
output  for  values  of  the  several  input  parameters  that  were  altered,  and 

figures  10  through  13  illustrate  results  for  different  values  of  TSOURC  and 

LIMIT. 

Changing  the  seed  (an  arbitrary  initial  input  value)  for  the  random  number 

generator  appeared  to  alter  the  overall  distribution  and  structure  of  the 

cloud  to  some  extent,  especially  the  height  of  the  highest  particle  (ZMAX). 

Table  6  indicates  a  spreading  of  the  vertical  distribution  only  to  slightly 
greater  heights,  except  fcr  the  top  few  particles.  The  maximum  number  of 
particles  per  layer  occurred  for  layers  6  and  7  (10  to  20  m),  the  same  as  the 
maximum  zone  of  table  2.  However,  ZMAX  differed  by  about  31  m,  and  the  over¬ 
all  structure  of  the  cloud  shown  by  the  statistics  T,  T,  ox,  was  altered. 

Two  other  runs  (not  shown)  with  seeds  of  5689  and  987654123  produced  much 
smaller  changes  in  the  overall  distribution  and  structure;  for  example,  the 
values  of  X,  Y  for  each  layer  more  closely  matched  those  of  table  2,  and  the 
ZMAX  only  differed  by  1  or  2  m.  However,  in  the  run  using  a  seed  of  987654123 
the  maximum  zone  was  about  10  m  higher.  These  results  suggest  that  when 
making  comparisons,  the  same  seed  should  be  used  to  avoid  confusion,  since  a 
change  of  the  seed  may  exaggerate  or  suppress  differences  produced  by  changes 
of  other  In put  parameters. 

Altering  the  cooling  rate  parameter  (COOL)  led  to  more  significant  changes 
(see  tables  7  and  8).  For  lower  values  of  COOL,  more  particles  rise  above  the 
zone  of  maximum  concentration,  even  though  most  remain  below  about  70  m  (below 
layer  18).  The  maximum  height  is  two  and  one-half  times  greater  for 
COOL  =  0.2  than  for  COOL  =  0.8.  An  equal  change  in  the  “drag"  coefficient 
(KFACT)  produced  more  dramatic  changes.  If  KFACT  =  0.8  (table  9),  the  verti¬ 
cal  profile  is  only  slightly  expanded,  and  ZMAX  is  only  about  15  m  higher  when 
compared  with  table  3  (instantaneous  release  without  buoyancy).  On  the  other 
hand,  if  KFACT  =  0.2  (table  10),  the  maximum  concentration  appears  to  lie 
between  about  45  and  80  m,  with  the  particles  strung  out  more  evenly  through¬ 
out.  Also,  a  zone  of  low  concentration  occurs  between  about  15  and  40  m,  and 
many  more  particles  rise  above  100  m.  The  listed  statistics  (  T,  T,  ox,  <0 
for  both  sets  of  tables  show  that  the  more  buoyant  clouds  (tables  8  and  10) 
spread  out  more  in  the  along-wlnd  and  crosswind  directions.  This  expansion  of 
the  cloud  probably  occurred  as  a  consequence  of  larger  scale  turbulence  and 
higher  windspeeds  as  height  increased.  Overall,  tables  7  through  10  suggest 
that  a  change  in  the  "drag"  coefficient  has  a  greater  effect  on  the  cloud  of 
particles  than  an  equal  change  in  the  cooling  parameter. 

If  we  change  the  proportion  of  particles  that  are  cooled  by  initial  contact 
with  the  ground  (or  are  considered  ingested  surface  material)  from  about  15 
percent  (table  2)  to  7.5  (table  11)  or  30  percent  (table  12),  the  main  effect 
seems  to  be  a  decrease  or  an  increase  in  the  number  of  particles  below  about 
10  m.  Otherwise,  the  vertical  profiles  remained  fairly  similar,  and  ^MAX 
changed  only  by  about  1  m. 
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Significant  changes  in  the  cloud  also  took  place  when  the  temperature  of  the 
source  (TSOURC)  was  varied  by  +30  K  (compare  tables  2,  13,  and  14).  The 
changes  resembled  those  that  appeared  when  KFACT  and  COOL  were  varied  (tables 
7  through  10).  "he  vertical  distribution  of  table  13  (TSOURC  =  383  K)  is 
similar  to  that  of  table  10  (KFACT  =  0.2)  in  the  spreading  of  the  particles 
more  evenly  with  height  and  in  the  placement  of  the  zone  of  maximum  concentra¬ 
tion  near  60  m.  The  profiles  are  roughly  similar  for  table  14  (TSOURC  = 
328  K)  and  table  9  (KFACT  =  0.8)  as  well.  Other  statistics  (for  example,  J, 
ox)  show  the  greater  extent  and  variability  of  the  cloud  when  the  source 
temperature,  and  therefore  buoyancy,  is  increased. 

Figures  10  through  13  present  in  graphical  form  the  output  appearing  in  tables 
11  through  14  (LIMIT  =  0.075,  0.30  and  TSOURC  =  388,  328  K).  These  figures 
nay  be  compared  with  figures  4,  5,  8,  and  9,  which  show  output  for  an  instan¬ 
taneous  release,  with  and  without  buoyancy,  for  X,  Z  and  Y,  Z  views.  Figures 
10  and  11  illustrate  the  greater  concentration  of  particles  near  the  surface 
when  more  particles  lost  their  excess  temperature  on  release  (that  is,  when 
LIMIT  increased).  Nevertheless,  the  overall  size  of  the  cloud  did  not  change 
greatly,  despite  an  increase  of  nonbuoyant  particles  on  release  from  7.5  to  30 
percent.  The  volume  of  the  cloud  changed  significantly  when  the  source  tem¬ 
perature  was  varied  by  +30  K,  as  shown  in  figures  12  and  13.  In  fact,  the 
cloud  was  not  much  bigger  than  for  the  nonbuoyant  case  (figures  5  and  9)  when 
TSOURC  dropped  to  328  K.  For  TSOURC  =  388  K,  the  cloud  expanded  signifi¬ 
cantly,  especially  in  the  vertical,  where  the  highest  particle  had  risen 
beyond  the  confines  of  the  expanded  vertical  scale  of  0  to  200  m.  Almost  all 
particles  beyond  a  downwind  distance  of  about  315  m  had  risen  above  a  height 
of  8  or  10  m.  This  elevation  of  the  bulk  of  the  cloud  above  the  surface  would 
permit  a  nearly  clear  horizontal  view  for  ground-based  systems. 

This  brief  preliminary  analysis  suggests  that  the  most  sensitive  input  coeffi¬ 
cient  of  those  considered  is  KFACT,  followed  closely  by  COOL.  Variations  in 
LIMIT  primarily  affect  the  lowest  10  m.  However,  different  seeds  for  the 
random  number  generator  may  produce  effects  large  enough  to  hide  or  enhance 
changes  in  the  particle  cloud  caused  by  varying  other  parameters.  The  varia¬ 
tion  of  TSOURC  resulted  in  smaller  differences  in  the  cloud  than  the  change  in 
KFACT.  However,  a  proportional  change  in  TSOURC  equal  to  that  in  KFACT  (30 
percent  or  about  108  K)  would  have  produced  the  greatest  alteration  in  the 
size  and  structure  of  the  cloud.  To  correctly  model  a  source,  we  would  have 
to  adjust  the  several  coefficients  at  the  same  time  as  the  source  temperature, 
a  complex  analysis  beyond  the  scope  of  this  preliminary  effort. 
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A  preliminary  buoyancy  routine  for  the  particle  model  originally  developed  by 
Ohmstede  and  Stenmark*  2  allows  for  the  realistic  simulation  of  the  transport 
and  diffusion  of  buoyant  material.  Future  improvement  to  the  algorithms 
should  include  the  computation  of  the  several  coefficients  based  on  atmo¬ 
spheric  and  source  properties,  replacing  the  present  need  for  input  values. 
In  line  with  this  goal,  a  more  complete  sensitivity  analysis  should  be 
performed  to  determine  in  detail  what  changes  in  which  coefficients 
significantly  affect  the  evolution  of  a  buoyant  cloud.  It  may  turn  out  that 
precomputed  values  for  specific  types  of  sources,  in  one  or  more  "look-up" 
tables,  may  suffice  for  most,  if  not  all,  applications. 

Addition  of  buoyancy  to  the  particle  model  reveals  some  Interesting  properties 
applicable  to  US  Army  systems,  such  as  the  much  greater  coverage  by  even  a 
modestly  buoyant  cloud.  Perhaps  more  Important  Is  the  possibility  that  a 
major  portion  of  a  strongly  buoyant  cloud  of  smoke  or  other  obscurant  may  lift 
off  the  surface  entirely,  allowing  personnel  and  systems  on  the  ground  to  see 
and  be  seen.  Thus,  the  best  screen  may  be  produced  by  a  modestly  buoyant 
cloud,  permitting  a  greater  coverage  per  round,  but  avoiding  the  creation  of  a 
"clear"  layer  near  the  surface. 

In  any  case,  the  model  with  buoyancy  should  permit  the  realistic  modeling  of 
maty  types  of  smokes,  dusts,  other  obscurants,  or  chemical  agents.  This 
buoyancy  routine  Is  one  step  In  a  series  of  Improvements  for  the  particle 
model  (such  as  evaporation  of  deposited  agent  and  settling)  that  are  planned 
over  the  next  several  years. 


LW.  0.  Ohmstede  and  E.  B.  Stenmark,  1980,  "A  Model  for  Characterizing 
Transport  and  Diffusion  of  Air  Pollution  In  the  Battlefield  Environment," 
Proceedings  of  Second  Joint  Conference  on  Applications  of  Air  Pollution 
Meteorology  and  Second  Conference  on  industrial  Meteorology,  flmer  Met  Soc, 
Boston,  MA 


2W.  D.  Ohmstede  and  E.  B.  Stenmark,  1981,  "Parameterization  of  the  Dispersion 
of  Battlefield  Obscurants,"  Proceedings  of  the  Smoke/Obscurants  Symposium  V, 
Adel  phi,  MD 
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.  Simplified  two-dimensional  sketch  of  a  sample  "cloud"  showing  the 
buoyant  "zone"  (within  heavy  line)  and  the  division  of  the  cloud 
into  layers.  In  the  model,  layers  from  the  surface  to  10  m  have  a 
thickness  (aZ)  of  2  m,  while  layers  above  10  m  have  a  aZ  of  5  m. 
Mean  (X)  and  standard  deviation  (+o  )  in  X  are  shown. 
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Curve  showing  the  portion  [1  -  ke  *  (1  -  f)  =  p]  of  the  difference 
in  temperature  (AT)  between  a  particle  and  the  ambient  atmosphere 
that  Is  retained  (p  *  AT)  for  a  given  time  step  and  layer.  A 
particle  at  the  center  Of,  Y)  will  not  lose  any  heat  by  entrain¬ 
ment.  Nine- tenths  of  AT  is  lost  In  the  "fully"  entrained  region 
(for  example,  where  ke  =  0.9  and  Y  >  Y  +  o  )  . 
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DECREASE  PARTM1.E  TEMPERATIKE  T«  ACCtlENT 
EOR  AIM  AHATIC  EXPANSION  (COOLIV.I. 

TP  ART  (\P)  =  TP  ART  (NP)  —  IMMMM  £  (ZZ(NP)—  ZPREVOP)) 


* 

I  COMPETE  TEMPERATERE  DIFFERENCE  I 
i  IIELT=TPART(NP>  —  TAIR  I 


TEST  IP  PARTICLE  HI  OY  ANT 
IIEl.Ti  II.  I 


TEST  IK  =»■  KIRST  TIAIK  STEP  KOR 

II  PARTlr|.E(NIM>>  I  AMI  NP  ^  Mji 


[SELECT  OAESSIAN  RANDOM  NCMBERl 


is  random  ntmber  •<  inpi  t 

VALEK  (I.IMIT)? 


COMPETE  PROPORTION  OK  TEMPERATERE 

oikkerenck:  mist 

IIPl'U  —IIVIIT  AIK  ITPAHT(VM) — TSI  HKl 


COMPETE  NEW  PARTICLE  TEMPERATERE 
TP  ART  |NIH=TPART  (NP) — HECK 


COMPETE  VAI.EE  OK  LAYER  INDEX  k 


TEST  IK  PARTICLE  OETSIIIE  OF 
CENTRAL  ELLIPSE  KOR  LAYER  k 


TEST  IK  £  3  PARTICLES  I 
LAYER  k  (if  <T  SS  <»> 


COMPETE  NEW  PARTICLE  TEMPERATERE 
TP  ART  (NP)  =  TP  ART  (NP)  —  KNTKAC^f  DELT 


COMPETE:  (»)<J^  O,  ,  IX— XT.  iT— Y)' 


(h)  PROPORTION  OK  TEMPER ATERE 
IHEKKKENCE  REMAINING  (XYKACT) 

(f)  NEW  PARTICLE  TEMPERATERE 

TP  ART  (NP)=  TP  ART  (NP)  — (I  -XYKA(T)  * 
ENTKAC  m  IIEI.T 


Figure  3.  Flowchart  of  the  min  buoyancy  routine.  ENTFAC,  GFACT,  KFACT,  and 
COOL  are  coefficients  Input  by  the  name  list  SPDATA.  The  limiting 
factor  (LIMIT)  for  the  Ingestion  of  ground  material  and  the  buoy¬ 
ancy  flag  (IBUOY)  also  are  Input  through  the  name  list.  Default 
values  are  provided  in  the  program.  Blocks  with  asterisks  along¬ 
side  were  part  of  the  program  before  buoyancy  was  added.  Brief 
definitions  of  Input  variables  are  found  In  table  1. 
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Figure  4.  Output  from  the  particle  Model  with  buoyancy  for  an  Instantaneous 
release.  See  table  2  for  other  Initial  conditions.  The  X,  Z  view 
Is  shown. 


Figure  5.  Output  from  the  particle  model  without  buoyancy  for  an  Instanta¬ 
neous  release.  See  table  2  for  other  Initial  conditions.  The  X,  Z 
view  Is  shown. 


0  50  100  150  200  250  300  350  400  450  5G0  550  600 

X 

Figure  7.  Output  fro*  the  particle  Model  without  buoyancy  for  a  quasl- 
contlnuous  release.  See  table  2  for  other  Initial  conditions.  The 
X,  Z  view  Is  shown. 
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Figure  12.  Output  from  the  particle  Model  with  a  buoyant.  Instantaneous 
release  after  150  s.  Initial  conditions  were  the  saw  as  for 
table  2,  except  that  TSOURC  *  388  K  In  the  upper  graph  and  328  K 
In  the  lower  one.  Compare  this  figure  with  figures  4  and  5. 


Figure  13.  Output  from  the  particle  model  with  a  buoyant.  Instantaneous 
release  after  150  s.  The  Input  values  were  the  same  as  for  figure 
12,  except  for  the  Y,  Z  view.  Compare  this  figure  with  figures  8 
and  9. 


TABLE  1.  INPUT  PARAMETERS  AND  VARIABLES 


Variable 

Defi nltion 

Default  Value 

I  BUOY 

Flag  for  use  of  buoyancy  routine. 

0  =  no  buoyancy,  otherwise  yes. 

0 

TSOURC 

Source  temperature  in  kelvin,  that  is,  temper¬ 
ature  of  a  particle  when  released. 

288 

(358) 

TSURF 

Surface  temperature  In  kelvin. 

288 

LAPSE 

Lapse  rate  ( -dT/dZ )  in  kelvin  of  the  ambient 
atmosphere  for  Z  >  HITE. 

0.002 

ADLAPS 

Lapse  rate  in  kelvin  for  Z  <  HITE. 

0.0098 

HITE 

Height  in  meters  of  boundary  between  bottom  and 
top  layers  of  ambient  atmosphere. 

50. 

KFACT 

A  "drag"  coefficient  used  to  simulate  drag  on  the 
additional  acceleration  arising  from  buoyancy. 

0.0 

(0.5) 

GFACT 

Coefficient  giving  the  proportion  of  AT  lost 
on  contact  with  the  ground.  Also  used  for  simu¬ 
lation  of  Ingestion  of  nonbuoyant  material  at  the 
surface. 

1.0 

ENTFAC 

Entrainment  coefficient  giving  proportion  of  AT 
lost  as  a  result  of  "full"  entrainment  (outside 
the  central  region  of  the  cloud). 

0.9 

COOL 

Coefficient  determining  the  rate  of  exponential 
cooling,  not  including  entrainment.  Ingestion,  or 
contact  with  the  ground. 

0.2 

(0.5) 

LIMIT 

An  adjustable  parameter  giving  the  proportion  of 
particles  assumed  to  lose  heat  on  initial  contact 
with  the  ground  or  to  consist  of  ingested  material 
(proportion  of  particles  assumed  to  be  nonbuoyant 
on  release). 

.15 

*The  values  used  for  the  first  sample  run  (table  2)  are  shown  In  parentheses 
If  they  differ  from  the  default  values. 


TABLE  2.  OUTPUT  FR<M  THE  PARTICLE  HOOEL  WITH  BUOYANCY* 


Miaber  of  Particles  and  Seconds  -  500  150.00000 


Layer  and  Particles/Layer 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


34 

35 
34 
21 
3C 
54 
50 
30 
41 
18 
22 
26 
23 
18 
17 

8 

12 

6 

4 

8 

2 

3 

1 

1 

0 

1 

0 

0 

0 

1 


ZMAX  =  133.40399 


XBAR 

XSIG 

YBAR 

YSIG 

1 

226.12803 

42.807256 

1.7038317 

14.633185 

2 

271.59393 

60.986114 

1.1277608 

26.299388 

3 

280.00192 

64.535798 

1.9918521 

27.908941 

4 

300.18904 

75.465451 

-7.5869053 

38.212638 

5 

307.06327 

70.255499 

11.510468 

30.360852 

6 

335.51123 

67.095635 

-3.7746255 

32.493101 

7 

354.62260 

56.518667 

-6.2751591 

39.322205 

8 

367.75878 

64.568353 

4.6168453 

36.489753 

9 

394.15417 

60.911937 

-4.3249745 

42.000307 

10 

400.00241 

59.647925 

14.224563 

40.426382 

11 

398.99990 

60.707179 

17.610682 

47.545129 

12 

403.51014 

45.475460 

2.7584606 

48.999747 

13 

405.11242 

61.446808 

14.941782 

43.275398 

14 

424.20483 

49.801628 

-5.6891095 

37.710571 

15 

440.47490 

43.553241 

-2.2888620 

48.448833 

16 

424.99738 

64.926211 

13.485317 

30.414957 

17 

443.67565 

49.304250 

5.5270264 

32.754390 

18 

453.24629 

36.696857 

-21.643807 

31.662502 

19 

503.84800 

35.041828 

19.483190 

60.730539 

20 

437.72428 

46.857068 

-16.834477 

41.142911 

♦Five  hundred  particles  were  released  Instantaneously  and  the  output  printed  for  a  model  time  of  150 
s.  K  *  layer  number,  N  -  number  of  particles  In  each  layer,  ZMAX  =  height  of  highest  particles,  XBAR 
and  YBAR  are  the  mean  X  and  Y  positions  In  each  layer,  and  XSIG  and  YSIG  are  the  respective  standard 
deviations.  The  radiation  Index  was  0  (neutral  atmospherel ,  and  the  wlndspeed  Increased  by  a  power 

law  (u  Zp,  where  P  *  0.2). 


TAB  E  3.  OUTPUT  FROM  THE  PARTICLE  MOOa  WITHOUT  BUOYANCY* 
Nwfcer  of  Particles  and  Seconds  *  500  150.00000 


Layer  and  Particles/Layer 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


151 

90 

85 

57 

35 

47 

17 

8 

8 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


ZMAX  «  37.626780 


JC 

XBAR 

XSI6 

YBAR 

YSI6 

1 

207.88064 

39.907544 

-.23669814 

16.534172 

2 

230.68118 

41.091720 

-4.9657560 

20.560721 

3 

254.51769 

44.913856 

2.9305366 

20.870197 

4 

244.83993 

51.983192 

-1.0774224 

27.008596 

5 

269.63817 

39.282068 

-6.0558610 

20.181514 

6 

272.87398 

46.402422 

3.5143632 

28.184915 

7 

280.70849 

43.206140 

9. IU 18606 

31.857800 

8 

292.30516 

77.137263 

5.7934824 

43.436589 

9 

321.59892 

35.167504 

5.3148208 

33.170277 

♦Other  Initial  conditions  w ere  the  same  as  for  Table  2. 


TABU  *.  OLTPUT  FOR  A  QUASI -CONTINUOUS  BUOYANT  RaEASE* 
Nuaber  of  Particles  and  Seconds  =  500  150.00000 
Layer  and  Particles/Laver 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


66 

63 

43 

48 

20 

60 

40 

38 

36 

26 

19 

11 

8 

7 

4 

4 

2 

2 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


ZMAX  =  82.562564 


Jt 

XBAR 

XSIG 

YBAR 

YSIG 

1 

203.91289 

40.543240 

-2.9302797 

15.934371 

2 

232.67127 

58.788380 

-1.3340548 

21.363602 

3 

256.63818 

61.909195 

-1.0017261 

21.788335 

4 

256.35050 

52.330463 

-1.4320815 

28.595548 

5 

268.69769 

51.392845 

.42345021 

34.460595 

6 

293.19323 

53.075129 

.44264154 

32.931817 

7 

323.12839 

63.354770 

-11.125333 

29.089018 

8 

319.31132 

55.657214 

-1.9971253 

31.997339 

9 

324.03738 

57.479212 

1.7913352 

38.498520 

10 

353.28620 

67.463490 

9.1417046 

42.499573 

11 

375.60770 

50.410950 

.80433124 

53.499415 

12 

369.79680 

38.312061 

7.7498828 

45.847095 

13 

350.25494 

42.086327 

-18.110951 

40.026799 

14 

407.66085 

41.494226 

4.8410304 

38.345977 

15 

372.06316 

26.390497 

-li. 392940 

15.433369 

16 

418.74789 

39.715394 

17.668413 

31.315602 

*The  release  had  a  burn  time  of  30  s.  Other  Initial  conditions  were  the  same  as  for  Table  2 


TABLE  5.  OUTPUT  FOR  A  qUASI -CONTINUOUS  NONBUOYANT  R&EASE* 
Number  of  Particles  and  Seconds  >  500  150.00000 
Layer  and  PartlclesAayer 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


156 

106 

80 

54 

36 

44 

12 

6 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


ZMAX  =  29.434073 


£ 

XBAR 

XSI6 

YBAR 

1 

180.08789 

37.176643 

-1.5898772 

2 

211.48266 

39.716743 

.80273122 

3 

213.01281 

43.708303 

.45429102 

4 

217.77686 

43.460036 

-2.0792077 

5 

227.74085 

47.664972 

-.65602818 

6 

245.40714 

39.353875 

3.8727165 

7 

233.07872 

41.601242 

.31544929 

8 

285.29615 

26.025213 

-1 1.072637 

9 

274.63539 

53.288725 

20.933601 

♦The  release  had  a  burn  time  of  30  s.  Other  initial  conditions  were  the  same  as 


YSI6 

13.864060 

19.676161 

20.982479 

21.588795 

22.076996 

28.872782 

30.466342 

24.251785 

34.780147 

for  Table  2 


TABLE  6.  OUTPUT  FOR  A  NEW  INPUT  SEED  FOR  THE  RANDOM  NUMBER  GENERATOR 
Nunfcer  of  Particles  and  Seconds  =  500  150.00000 


Layer  and  Partlcles/Layer 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


33 

26 

23 

27 
20 
62 
53 
38 

28 

32 
26 

33 
19 
21 
11 

8 

11 

9 

5 

3 

4 
2 

5 
1 
0 
0 
0 
0 
0 
0 


ZMAX  =  101.88895 


JC 

XBAR 

XSIG 

YBAR 

YSIG 

1 

216.64240 

38.954814 

-9.2285895 

14.157842 

2 

278.65093 

75.914277 

.37619312 

27.219905 

3 

272.92013 

52.755116 

-5.0258740 

22.424229 

4 

290.98133 

79.362779 

3.0338036 

26.095520 

5 

339.85266 

57.438826 

7.6432076 

36.281323 

6 

352.495't4 

68.215305 

-5.9814332 

41.872091 

7 

350.01913 

59.860967 

-. 245u4836 

43.039150 

8 

386.38350 

63.434380 

2.0324648 

43.623247 

9 

398.58367 

49.273471 

-3.3974730 

39.315310 

10 

381.82377 

55.185519 

-8.2368256 

38.107791 

11 

399.96571 

51.433865 

-13.947372 

39.464909 

12 

401.50166 

53.623233 

2.9903516 

39.504470 

13 

411.07472 

51.133930 

-11.804703 

47.162509 

14 

417.34556 

46.915832 

2.0568820 

41.426896 

15 

434.30298 

42.143060 

31.237015 

34.400537 

16 

401.00705 

44.517277 

3.5167478 

38.070114 

17 

451.03362 

38.554142 

11.626332 

48.211498 

18 

425.62295 

55.574398 

4.3184916 

64.363358 

19 

429.84867 

31.470472 

3.9468845 

33.277129 

21 

444.61465 

52.385838 

-2P. 622247 

60.586889 

23 

448.30622 

52.457870 

.85011664 

40.283696 

•Initial  conditions  were  the  same  as  for  Table  2,  except  that  the  seed  for  the  random  number  generator 
was  changed  to  123987654  from  321456789. 


TABLE  7.  OUTPUT  FOR  A  LARGER  INPUT  VALUE  OF  THE  COOLING  RATE  COEFFICIENT* 


N'Mtter  of  Particles  and  Seconds  -  500  150.00000 


Layer  and  Particles/Layer 
K  N 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


39 

40 

36 
28 

37 
67 
51 
47 
39 
27 
25 
22 
11 
14 

8 

4 

1 

1 

2 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 


ZMAX  =  96.481231 


JC 

MAR 

XSI6 

YBAR 

YSIG 

1 

235.06515 

50.353797 

-.15403188 

16.185288 

2 

268.95041 

56.422684 

-1.2455061 

28.113550 

3 

284.73095 

66.586268 

8.5215377 

22.962459 

4 

314.15090 

74.422506 

-7.8220468 

39.816657 

5 

300.94077 

69.086676 

-4.2150003 

29.952330 

6 

341.01091 

63.098648 

-.75165023 

38.401177 

7 

352.20020 

58.219295 

.33320645-001 

35.462511 

8 

377.75522 

58.161145 

8.0532399 

40.868855 

9 

383.34976 

60.707276 

-5.3917148 

40.356915 

10 

383.98060 

48.091472 

1.6648052 

34.413349 

11 

399.66768 

45.816252 

9.0257194 

43.331216 

12 

393.15466 

52.714852 

7.0382617 

46.285177 

13 

398.81858 

60.079050 

-3.5034595 

41.297105 

14 

418.17549 

53.758921 

-.57404654 

42.745411 

15 

408.18529 

71.602293 

25.655352 

36.193396 

16 

443.80283 

70.787038 

1.4306523 

14.406935 

♦Initial  conditions  were  the  same  as  for  Table  2,  except  that  the  coefficient  for  cooling  rate  COOL 
0.8  Instead  of  0.5. 


TABLE  8.  OUTPUT  FOR  A  SMALLER  INPUT  VALUE  OF  THE  COOLING  RATE  COEFFICIENT 


Ntafcer  of  Particles  and  Seconds  =  500  150.00000 


Layer  and  Particles/layer 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


34 

22 

26 

26 

23 

51 

42 

32 

30 

28 

14 
18 
18 

8 

15 
7 

14 

4 

6 

4 

5 
7 
3 
7 

10 

3 

2 

5 

2 

2 


ZMAX  =  238.79427 


JC 

XBAR 

XSI6 

YBAR 

YSIG 

1 

230.51976 

47.944644 

1.6986663 

14.866672 

2 

267.12107 

66.148516 

3.6396527 

26.334696 

3 

273.09214 

59.478798 

-.16151821 

28.434213 

4 

298.73883 

77.748878 

-.82682638 

23.112719 

5 

308.60403 

79.734155 

13.205221 

39.886557 

6 

336.81034 

67.851704 

-8.8032271 

34.891039 

7 

354.77332 

60.746236 

1.3272103 

43.218959 

8 

364.01276 

59.910972 

-.56122220-003 

34.024921 

9 

393.46827 

66.824300 

10.292469 

36. 723G27 

10 

391.68689 

53.430012 

-6.8423545 

39.178792 

11 

392.48566 

66.398618 

25.550022 

41.634145 

12 

402.08502 

31.295672 

5.5655505 

46.938416 

13 

404.04905 

67.394634 

10.002160 

54.289927 

14 

427.83150 

67.799402 

-15.937131 

29.452288 

15 

450.36797 

49.994754 

-9.4677540 

45.154819 

16 

414.58881 

70.497709 

8.5421247 

28.404871 

17 

449.16714 

45.570792 

7.7554339 

55.467068 

18 

462.23301 

82.339000 

-13.348862 

32.210813 

19 

413.54.89 

56.321607 

1.6880149 

40.532016 

20 

458.26082 

47.841753 

-.16808838 

21.588892 

21 

468.63512 

40.522322 

-13.31C712 

44.355479 

22 

460.78326 

39.562319 

34.614296 

31.452911 

24 

487.63681 

37.618320 

-33.881571 

43.015985 

25 

496.75961 

37.598390 

25.702019 

49.866966 

28 

472.36997 

39.783840 

-1.3407461 

41.217841 

•Initial  conditions  were  the  same  as  for  Table  2,  except  that  the  coefficient  for  cooling  rate  COOL 
0.2  Instead  of  0.5. 


TABLE  9.  OUTPUT  FOR  A  LARGER  INPUT  VALUE  OF  THE  ‘DRAG*  COEFFICIENT’ 
Nuafcer  of  Particles  and  Seconds  =  500  150.00000 
Layer  and  Particles/Layer 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


122 

101 

72 

54 

39 

50 

27 

15 

13 

3 

1 

2 

0 

1 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 


ZMAX  =  52.193466 


Jt 

XBAR 

XSIG 

YBAR 

YSIG 

1 

219.60612 

45.651437 

-1.8275875 

17.331466 

2 

249.48536 

48.307285 

-1.7782056 

19.623469 

3 

259.43560 

55.152477 

.68552608 

19.343807 

4 

263.34941 

55.815402 

-4.7818034 

23.572031 

5 

275.85178 

52.030853 

13.226504 

32.250138 

6 

299.96722 

50.792821 

.73963130 

31.605940 

7 

313.35028 

61.735186 

5.0209183 

34.418818 

8 

321.30085 

51.416823 

8.3418875 

32.252311 

9 

327.71467 

52.172241 

-2.8729145 

25.933355 

♦Initial  conditions  vie  re  tn?  same  as  for  Table  2,  except  that  the  "drag"  coefficient  KFACT 
Instead  of  0.5. 


TABLE  10.  OUTPUT  FOR  A  SMALLER  INPUT  VALUE  OF  THE  'DRAG*  COEFFICIENT 
TMfcer  of  Particles  and  Seconds  *  500  150.000C0 
Layer  and  ParticlesAayer 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


23 

16 

11 

12 

10 

11 

7 
4 
6 

8 
9 

10 

23 
22 

24 
29 
18 
24 
28 

15 

16 
8 

18 

15 

15 

11 

7 

8 
7 
9 


ZMAX  «  249.75252 


t 

XBA8 

XS1G 

YBAR 

TSIG 

1 

207.23906 

30.108492 

.76983353 

13.851210 

2 

226.53797 

40.575497 

5.8508274 

23.383023 

3 

224.92763 

38.887639 

9.8974869 

21.878350 

4 

268.45663 

76.682816 

-5.0382897 

21.723365 

5 

265.76254 

71.270983 

20.675223 

36.214669 

6 

320.55933 

104.24301 

-19.324917 

31.998740 

7 

316.26085 

68.108962 

8.9293455 

60.077665 

8 

413.80284 

44.585709 

7.1943512 

19.716362 

9 

430.53845 

77.969886 

-25.170937 

56.882793 

10 

403.93951 

65.835682 

4.0253773 

17.868464 

11 

430.26085 

46.449254 

11.066533 

23.462560 

12 

401.17316 

67.593518 

-15.992053 

49.093642 

13 

434.09926 

69.394693 

.31037183 

40.856453 

14 

433.12608 

60.789670 

1.0315800 

55.021161 

15 

442.89719 

69.289131 

.24497608-002 

52.081959 

16 

442.17851 

53.886827 

11.040929 

36.518641 

17 

438.40559 

54.007897 

-7.7791709 

40.380655 

18 

459.40825 

65.158377 

-3.1647065 

46.580976 

19 

444.17250 

65.736972 

16.175095 

45.134769 

20 

476.40837 

59.365977 

-10.691096 

55.750010 

21 

478.36372 

68.941476 

18.951031 

53.055208 

22 

474.54899 

45.463903 

4.7218918 

59.035519 

23 

458.10590 

46.725339 

-.71359132 

56.178855 

24 

463.94514 

57.187622 

16.368306 

49.226245 

25 

479.98733 

58.842556 

-2.4120045 

46.635083 

26 

504.94154 

54.540610 

-17.470158 

30.307095 

27 

535.87697 

44.879323 

16.627070 

52.780869 

28 

470.89333 

59.117977 

-.73022604-001 

47.619286 

29 

517.47522 

45.468882 

-11.392206 

31.984057 

30 

488.31564 

62.752926 

1.5031427 

45.483222 

•Initial  conditions  were  the  same  as  for  Table  2,  except  that  the  "drag"  coefficient  KFACT 
instead  of  0.5. 
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ZMAX 


TABLE  11.  OUTPUT  FUR  A  SMALL  ER  INPUT  VALUE  OF  THF  PROPORTIONALITY 

COEFFICIENT  FOR  ENERGY  LOSS  ON  RELEASE* 


Nurrfcer  of  Particles  and  Seconds  =  500  150.00000 


Layer  and  Parttcles/Layer 
K  N 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


27 

29 
31 
13 

28 
53 
51 
34 
43 
21 
27 

30 
23 
2? 
15 

9 

12 

6 

5 

10 

2 

3 

2 

0 

1 

1 

1 

0 

0 

1 


132.41249 


£ 

XBAR 

KSIG 

YBAR 

YSIG 

1 

232.20476 

43.430585 

1.9993522 

15.105388 

2 

277.69165 

60.444198 

.55132417 

28.430081 

3 

295.27971 

60.136797 

1.0576700 

27.804968 

4 

327.36686 

65.699719 

-7.5013586 

45.276865 

5 

324.55400 

60.164117 

11.025974 

29.213091 

6 

339.80689 

64.493114 

-2.7007304 

32.387435 

7 

359.74087 

51.983171 

-3.4740485 

39.016682 

8 

373.37694 

65.583598 

1.4700783 

37.026200 

9 

389.60714 

61.239187 

-3.7993517 

41.750751 

10 

398. 1 5 1 1C 

62.083676 

2i. 165695 

48.823731 

11 

399.96888 

62.165528 

7.3675197 

45.008676 

12 

397.79013 

55.881127 

3.9715701 

51.875017 

13 

409.09315 

56.759523 

1C. 925770 

42.744894 

14 

429.09528 

51.719596 

-1.81202.9 

38.544679 

15 

437.16161 

45.423700 

-6.5732731 

48.682740 

16 

420.67117 

62.405585 

5.1215112 

37.910660 

17 

443.54973 

49.495207 

5.5307787 

32.766469 

18 

439.21038 

37.008867 

-19. 599832 

29.870799 

19 

500.73920 

31.242749 

5.9045761 

60.749024 

20 

423.74015 

53.842201 

-17.569979 

38.361137 

♦Initial  conditions  were  the  same  as  for  Table  2,  except  that  the  :oefficient  (LIMIT)  to  determine 
the  proportion  of  particles  that  lose  their  buoyant  energy  on  release  =  0.075  instead  of  0.15. 


ZMAX  = 


TABLE  12.  OUTPUT  FOR  A  LARGER  INPUT  VALUE  OF  THE  PROPORTIONALITY 

COEFFICIENT  FOR  ENERGY  LOSS  ON  RELEASE* 


Nunfcer  of  Particles  and  Seconds  =  500  150.00000 


Layer  and  Partlcles/Layer 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


49 

49 

41 

33 

28 

52 

44 

29 

30 
16 
20 
24 
17 
12 
15 

9 

8 

6 

7 

4 

1 

2 

1 

1 

0 

1 

0 

0 

0 

1 


134.50200 


K 

XBAR 

XSIG 

YBAR 

YSIG 

1 

215.29210 

38.121687 

2.0198779 

14.959438 

2 

248.79433 

56.792517 

-1.1165867 

24.652422 

3 

263.77854 

63.979361 

4.4112783 

24.934906 

4 

284.14340 

72.17205E 

-1.7914171 

32.928373 

5 

295.13359 

69.025030 

11.412329 

27.633392 

6 

329.00245 

67.538606 

-3.5798484 

23.806411 

7 

352.07703 

56.583353 

-8.3843718 

38.219877 

8 

370.14833 

70.038039 

-.78563196 

35.489929 

9 

382.18621 

61.3,0192 

-9.0408024 

42.093504 

10 

399.44917 

59.435907 

22.699064 

31.429122 

11 

406.71499 

62.883991 

14.277443 

44.355289 

12 

403.36322 

38.965425 

8.9888580 

41.208457 

13 

406.17208 

63.145325 

13.363133 

45,513962 

14 

425.90897 

46.887437 

-21.211224 

30.851888 

15 

438.39655 

43.893184 

-6.4368755 

47.591279 

16 

437.92161 

45.602263 

25.190234 

31.226167 

17 

432.65161 

51.762464 

1.3197374 

39.186940 

18 

468.77160 

24.121178 

-21.138547 

31.337169 

19 

448.78312 

64.933961 

-4.1034228 

57.333452 

20 

464.00409 

43.675293 

-18.487511 

38.374288 

•Initial  conditions  were  the  «ame  as  for  Table  2,  except  that  coefficient  LIMIT  =  0.30  instead 


TABLE  13.  OUTPUT  FOR  A  MARKER  SJURCE  TEMPERATURE* 


Ntaber  of  Pai  .cles  and  Seconds  =  500  150.00000 


Layer  and  PartlcTes/Layer 
K  N 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


24 

16 

13 
15 
11 
12 
12 
10 
15 

27 

24 

25 

28 
24 
31 
22 
15 
12 
20 
18 
10 

14 
8 

17 

7 

5 

6 

7 
6 

8 


ZMAX  =  229.37799 


£ 

WAR 

XSIG 

TBAR 

TS16 

1 

213.85997 

43.807486 

3.0373522 

17.518930 

2 

226.53788 

40.576612 

5.8508322 

23.383015 

3 

248.50162 

83.282186 

4.1525235 

31.691045 

4 

291.24471 

84.524636 

.39829416 

27.907353 

5 

272.18805 

71.228308 

16.292689 

32.120471 

6 

312.79946 

69.875122 

-9.6100578 

20.545622 

7 

359.23529 

96.711276 

-17.486562 

51.923285 

8 

399.03793 

41.688203 

5.7962136 

14.417962 

9 

393.86873 

51.552547 

-11.061272 

46.576862 

10 

393.96765 

72.508271 

-2.1403582 

44.631855 

11 

412.66687 

63.010933 

13.269316 

41.804791 

12 

425.66868 

61.094083 

-8.1751883 

47.497755 

13 

425.93071 

54.941509 

6.2261650 

37.787997 

14 

417.80804 

54.287780 

2.0944006 

44.881759 

15 

430.10677 

72.481894 

6.1047679 

46.239466 

16 

456.54841 

55.673527 

-.13351406 

47.159827 

17 

442.55015 

69.271496 

8.1654725 

50.735657 

18 

476.55069 

53.501912 

8.9004390 

53.402852 

19 

441.78776 

33.630832 

8.6865519 

54.029192 

20 

450.93009 

52.831656 

5.9436646 

54.199925 

21 

462.43130 

85.652601 

2.1919701 

37.677923 

22 

478.51771 

50.855730 

-19.476986 

37.613358 

23 

504.90457 

41.467500 

2.1901492 

60.510386 

24 

468.45487 

58.499299 

2.4698740 

34.369756 

25 

500.46824 

38.767469 

16.701744 

59.297643 

26 

460.06360 

58.525468 

9.3525681 

11.556052 

27 

489.43712 

40.854314 

12.161410 

52.541036 

28 

492.37316 

44.472487 

7.4405220 

42.754703 

29 

479.40260 

36.895714 

17.877722 

60.877487 

30 

517.40580 

30.527578 

-2.1232033 

43,955000 

•Initial  conditions  were  the  same  as  for  Table  2,  except  that  the  temperature  (TSOURC)  of  the  source 
388  K  Instead  of  358  K. 
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TABLE  14.  OUTPUT  FOR  A  COOLER  SOURCE  TEMPERATURE* 
Nuntoer  of  Particles  and  Seconds  =  500  150.00000 
Layer  and  Particles/layer 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


106 

70 

68 

48 

41 

65 

37 

24 

15 

15 

6 

0 

2 

2 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


ZMAX  =  63.257215 


JC 

XBAR 

XSIG 

YBAR 

YS1G 

1 

232.18460 

51.233685 

-3.0361313 

19.65823! 

2 

258.37076 

54.794288 

-.35588692 

21.333701 

3 

268.98148 

55.993895 

3.7329751 

24.358226 

4 

284.62764 

60.023509 

-6.2069080 

22.484486 

5 

294.72088 

64.830645 

2.2846254 

27.415254 

6 

314.67526 

52.254037 

4.5368927 

38.281246 

7 

323.97559 

59.865075 

1.5587415 

38.236937 

8 

346.47226 

51.149264 

-2.0893207 

35.203655 

9 

373.06096 

48.721718 

6.3391777 

32.910414 

10 

349.62811 

63.29370. 

3.1034001 

38.076659 

11 

371.64274 

52.428091 

5.3594352 

13.314982 

♦Initial  conditions  were  the  same  as  for  Table  2,  except  that  TSOURC  =328  K  instead  of  358  K 
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